Distribution function of persistent current 
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We introduce a variant of the replica trick within the nonlinear sigma model that allows calculating 
the distribution function of the persistent current. In the diffusive regime, a Gaussian distribution 
is derived. This result holds in the presence of local interactions as well. Breakdown of the Gaussian 
statistics is predicted for the tails of the distribution function at large deviations. 
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Introduction. - A striking manifestation of quantum 
mechanics in the mesoscopic physics of electrons is 
that an equilibrium persistent current (PC) can flow in 
normal-metallic rings threaded by a magnetic flux 
In the diffusive regime, this property arises from a flux- 
dependent interference contribution to the electron den- 
sity of states Q that survives in presence of the impurity- 
induced static potential disorder. The amplitude of PC 
is quite small and varies strongly from sample to sample: 
The theory for non-interacting electrons [3, 0] predicts 
that the ensemble average vanishes, while the typical am- 
plitude is $o-periodic ($o = h/e is the flux quantum) 
and scales as 7t yp ~ e/773, where td is the diffusion time 
along the ring. Thus, 7t yp ~ In A in micrometer-size rings 
made of conventional metals. Electron-electron interac- 
tions were predicted to induce a $o/2-periodic average 
current of order I av ~ X ee/rr), where A e ff is an effec- 
tive coupling constant In superconducting rings, this 
would yield a diamagnetic average current due to super- 
conducting fluctuations well above the superconducting 
critical temperature @. 

Early experiments measured the PC by detecting the 
small magnetic field that they produce with supercon- 
ducting quantum interference devices (SQUIDs). Both 
the $n/2-periodic average current in large ensembles of 
rings [3, [8| and $o-P Gr iodic current in single rings and 
small ensembles of rings were recorded; a low-flux 

diamagnetic response was observed in 0, [ll|. However, 
it was not always possible to reconcile their amplitude 
with the theories described above. This could be due to 
subtle effects related to canonical vs. grand canonical av- 
eraging in isolated rings [HI], magnetic impurities [Hj], 
or the great sensitivity of the PC to its electromag- 
netic environment [Hj]. A recent experiment on several 
rings addressed separately with a scanning SQUID micro- 
scope showed however a good agreement with the non- 
interacting theory for the typical current and no sign of 
average current 



An experimental breakthrough was made recently by 
measuring PCs with high-precision cantilever torque 
magnetometry [jjj ■ Notably, the technique allows better 
sensitivity and works under large magnetic fields (com- 
pared with the SQUID technique). The PC could thus be 



recorded over a huge number of flux periods. The ampli- 
tude of the typical PC was found to be in good agreement 
with the prediction for non-interacting electrons and no 
average current was detected. 

The aim of the present work is to address the distribu- 
tion function of PC that seems within reach of the new 
experimental technique. We demonstrate that in the dif- 
fusive regime the statistics is Gaussian. It justifies char- 
acterizing the PC with its first two cumulants only. The 
Gaussian distribution carries on in the presence of local 
interactions as well. To derive the result, we make use of 
a replica trick in the nonlinear sigmal model that allows 
obtaining the distribution function at once. It provides 
an alternative to the evaluation of all cumulants order by 
order [lit [lij . It could serve as a useful starting point to 
address questions such as canonical vs. grand canonical 
averaging, localization and cross-over to non-Gaussian 
statistics. 

Gaussian distribution for the persistent current of non- 
interacting electrons. - The PC flowing in a metallic ring 
pierced by a magnetic flux <I>, 
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OF 



(1) 



is related to the free energy F($) = — fcTln Z$, where Z$ 
is the partition function. The PC distribution function 
at given $ is the probability density for /(<&) to be equal 
to I, 



P(l) = {6 (1 -/(*))>, 



(2) 



where the brackets denote an ensemble average over dif- 
ferent disorder configurations in the ring. Using the iden- 

and the definition of the dcriva- 



tity 6(s) = J (dx/2Tr) 
tive, we express Eq. as 



P(l) = lira / — e 



iXx 



ryn ryn \ 



(3) 



where n = —n' = —ixkT/(<& - $'). 

The disorder-averaging for free fermions in Eq. (|3l) can 
be performed with a variant of the replica trick [201 ] in 
the fermionic non-linear sigma model 2l|,l22|. For this, 
we consider that the system is formed of n replicas with 



2 



flux $ and n' replicas with flux $' (n and v! integers) 
and wc evaluate 



(ZlZl.) = \ VQe 



-S[Q] 



(4) 



Here, Q is a matrix field acting in the direct product of 
the replica space of dimension n + n' , infinite Matsub- 
ara energy space, two-dimensional Gorkov-Nambu space 
(Pauli matrices r a ), and two-dimensional spin space 
(Pauli matrices <7fc). The Q matrix obeys the nonlinear 
constraint Q 2 = 1 and the charge conjugation symme- 
try Q = Q = T\aiQ T 02^1 ■, where Q T stands for the full 
matrix transposition. The action of the model is [22, 23 1 



normal to the ring, is used. We introduce a coordinate 
y along the ring measured in units of L, a flux vector 
4> = {0a}a in replica space measured in units of $0j with 
components (j) a = <t> = $/*o (<t>a =(/>' = $'/*o) f° r 
1 < a < n (n < a < n + n'), and s = e/E c , where 
E c = hD/R 2 is the Thouless energy related to diffusion 
time (E c = Air 2 h/T D ). Then, Eq. © reads 



S[Q] 



J dyTr [{d y Qf - 167T 2 er 3 Q] 



(6) 



S[Q] = yJ drTl ' i D ^ 2 - 4er 3^] 



(5) 



where S is the mean level spacing in the ring and d y = 

fie 

0) = Q(y = 1) and 



Here, v is the single-particle density of states per spin, 
D is the diffusion coefficient in the metal, d = V + 
(ie/H)[AT3, .] includes the effect of a vector potential as- 
sociated with the magnetic field B = rotA (different in 
each replica) , e is a fermionic Matsubara energy, and the 
trace 'Tr' is taken over all spaces of the Q matrix. 

Let us consider a quasi one-dimensional circular ring 
with length L = 2irR (R is the radius of the ring). The 
circular gauge A = where it is a unitary vector 



J 



V y + 2wi[(f)T3, .]. The single- valuedness of the Q field 
fixes the boundary conditions: Q{y 
V y Q(y = 0) = V y Q(y = l). 

In the metallic regime, the ring's conductance mea- 
sured in units of the conductance quantum, g = 
E c /(2w5), is large. Thus, we can evaluate Eq. (j4} within 
the saddle-point approximation. The Qq field that mini- 
mizes the action (|6]) is proportional to unity in spin and 
replica spaces, and diagonal in Matsubara space with 
value Qo( e ) = T 3 s ig n (£)- The action at the saddle point 
is (n + n')So where 5 = —(2irE c /5) J2 E k|; it does not 
depend on the flux. 



In order to study fluctuations near this saddle point we write matrices close to Qo as 

Q = Q {l + W + W 2 /2 + ...). 



(7) 



The constraints on the Q field imply that {Qo, W} = and W = —W; the requirement of convergency of the sigma 
model on the perturbative level implies that W' — —W. Then, we decompose the W field into its elements in Nambu 
and spin spaces, and its Fourier components: 



ea.e'b 



(y) 




ea.e'b ea,e'b 

-w k ,{~ p) * s k w k{p ,i 

e'b.ea ^ K —e'b,—ea 



1 0ke 



2iripy 



(8) 



Nambu 



where a- = ± for k = 1,2,3 and k = 0, respectively, and the components W k ^, b (W k ^, b ) only exist at ee' < 



(ee' > 0). An independent set of complex integration variables is then obtained with W k ^, b and W^ P) , b at a > b, 

and W*$ a and W^ a at e > 0. 

Expanding the action (|6|) near the saddle point up to quadratic terms in W, one finds: 

s (2) [w] = ^£EE£{[N + + (p+^ - « 2 ] \ w s%>»\ 2 + [N + k'l + (p + ^ + 4> b ) 2 ] \w k ^, b \ 2 } (9) 

ee' ab k p 



Mp) 



The Gaussian integration over the W field is then 
straightforward and yields: 

(Z n Z n ' r ) = e ~ n2s ^^~ n ' 2E ^' ^')-^in'E(4>,4>') /^q\ 

where we have omitted a factor which is equal to 1 in the 



replica limit n' = —n, and 
~(«M')=4£ £ £ln[ £ + e ' + (p + 0- S 0') 2 ]- 

p e,e'>0s=± 

(11) 

The replica trick now consists in assuming that 
Eq. (|10l) can be analytically continued to pure imaginary 



3 



variables n and n' , with n' 
([TO]) into ©, wc find 



-n. Then, inserting Eq. 



P(l) 



dx 
2^< 



iXx , 



-I 2 /(2/ t 2 vp ) 



/ \/ 27rI typ' 

(12) 

where 4 p = -2(fcT/$ ) 2 d 2 S(</>, ^')/W)lf=f That 
is, the distribution function of the PC is Gaussian with a 
zero mean value and typical value I typ . Using the Poisson 
summation rule, we can convert the sum over p in S into 
an integral and finally obtain, at zero temperature, 



24£ c 2 



^4 sin2 (27rg$/$o). 

,- i y 



(13) 



Thus, the PC has typical amplitude E c /§q ~ e/To. The 
distribution function (|12p is in agreement with the known 
results [4J for the average current (/($)) and its cumu- 
lant ((/($)/($')» = (/(*)/(*')> " ( / ( $ ))( / (*')> in non- 



interacting diffusive rings, with I typ = ((iX^) 2 }} 1 / 2 [24 1. 

This section contains the main result of this article, 
Eq. (fT2l) . In the following we illustrate several directions 
where it can be extended. 

Distribution function for the harmonics. - It may be 
more convenient experimentally to characterize the flux- 
current relation by its harmonic content. In this section, 
we show that the distribution function for the harmonics 
in the diffusive regime is also Gaussian. 

Due to time-reversal symmetry and flux-periodicity, 
the current-flux relation, 



J($) = ^/ g sin27r#, 

9=1 



(14) 



is fully characterized by its harmonics I q . The distribu- 
tion function for the harmonics, 



P q (J) = (S(J-I q )), 



(15) 



can also be determined with a replica trick. Indeed, using 
Eqs. (dJ, (fT4"[) . and integrating by parts, we first note that 



1/2 

(87rg/$ ) Jo # cos(27r -q<f))F(<fi). By definition of 



the integration, it also reads: 



JV 



I q = lim (47r?/JV$ ) Vcos(27r?^)F(^), (16) 



JV 



e=i 



where <p£ — £/ (2N). Now, inserting the representation of 
the delta- function and Eq. (fTB"]) into ([T5|. we find 



dx 

P q {J) = lim / — e 



N 



ijx 



(17) 



1=1 



where ni = iAirxqkT cos(27r<7</^) / '($>qN). 

The average over the disorder can also be performed 
within the fermionic sigma model by considering that the 



system is formed of ng replicas (ri£ integer) with flux 4>£ 
(1 < £ < N). In the saddle point approximation, one 
would find as a generalization of Eq. (pi 



exp 



E 



(18) 



Taking the replica limit, one again obtains a Gaussian 
distribution P q (J) ~ cxp(— J' 2 /2((/ 2 ))) for the harmon- 
ics, with zero average value and variance 



32fc 2 T 2 



1/2 



d(f>d(f> sin(27rg0) sin(27r<7<//) 



(19) 



In particular, ((J 2 )) = 96 / (n 2 $lq 3 ) at T = 0, in agree- 
ment with Eq. ([TBI) . 

Interactions. - The effect of electron-electron interac- 
tions can also be taken into account. To be specific, we 
consider the case of attractive, local pairing between elec- 
trons with opposite spins that was theoretically debated 
after the early experiments on PC. Then, the action ([5]) 
should be supplemented with an interaction term 22j, 



E 



16 

a 

+ (trr 2 (3 Ta , Ta ) 2 ]] , 



J drdr [(trTiQ ra , ra ) 5 



(20) 



where A is the Bardcen-Cooper-Schricffcr coupling con- 
stant, t is imaginary time, and the trace 'tr' is taken over 
spin and Nambu spaces. 

Above the superconducting critical temperature, the 
action can be evaluated in the Gaussian approxima- 
tion near the metallic saddle point Qq. For the one- 
dimensional ring, Eq. (|20[) results in an interacting con- 
tribution adding to ©: 

sPhw] = - 4Xn2T v V w 0(p \ w°, ip) : 

p,u>,a £,e'>0 

(21) 

where u is a bosonic Matsubara energy (also measured 
in units of E c ). Gaussian integration over the W field 
including Eqs. ([5]), ([2"T]) can be performed; it yields 

(Z$Z$,) X = {Zlzi) x = Q e~ n ~"^- n '—^'\ (22) 



where 



Sint(0) =E ln 



i-E 



{AXnkT/E c ) 
2e - \u\ + (p + 20) s 



(23) 



Here, fl fixes the energy bandwidth around the Fermi 
level over which pairing is effective. By introducing the 



4 



critical temperature T c ~ (1.14Q/k)e 1 ^ A , one gets 



, T ,/l U + (p + 2(j)) 2 



-r 



.(24) 



where ip is the digamma function. Inserting Eq. ([22]) 
into ([3]), one again finds that the distribution func- 
tion for PC is Gaussian, with average value (/($)) = 
— (fcT/$o)9S; n t/90 and the same variance as in the non- 
interacting case. The average current was discussed 
in Ref. |6(, it is <£>o/2-periodic with amplitude / av ~ 
Aeffe/rD where A e fi ~ \n~ 1 {E c /kT c ) at T c < T <C E c . 

Discussion. - We first note that spin and orbital effects, 
such as the penetration of the magnetic field within rings 
with finite thickness [25| are important for a quantita- 
tive comparison with the experiment 17J . Taking these 
effects into account within our formalism can be done 
easily; it would not change the prediction of a Gaussian 
distribution in the diffusive regime. 

On the other hand, the Gaussian statistics clearly fails 
in the insulating regime, at g < 1. Actually, its alteration 
is expected already at large, but finite g, in relation with 
the Anderson localization phenomenon. A similar ques- 
tion on the statistics of the - dissipative - conductance 
of diffusive wires was addressed |26[. Log- normal tails 
in the probability distribution were predicted at large 
deviations from the average conductance. However, the 
present case differs by the fact that PC is a thermody- 
namic quantity. 

To estimate the range of validity of the Gaussian statis- 
tics for PC, we expand the action §5§ in vicinity of the 
metallic saddle point Qq up to fourth order terms in the 
field W. Then, we evaluate the generated terms pertur- 
batively with the Gaussian action. As a result, we found 
that the leading correction to the integrand in Eq. (fT2")) 
arises in order oc x 3 I^ yp /g, consistent with the recently 
derived third order cumulant (28[. The same way, we 
also obtain that n-th order cumulants scale as I™ yp /g n ~ 2 
at n > 3. Subsequently, this implies that the Gaus- 
sian distribution is not reliable at large deviations, when 
l^-l <L ff 1 ^ 3 ^typ 3> ityp- A more detailed investigation of 
the behaviour of P(X) at large deviations is left for future 
study. 

In the absence of interactions, the average current 
vanishes. However, this result is an artefact of the 
grand canonical averaging tacitly performed here. When 
(canonical) averaging is done with keeping the number of 
electrons constant in the ring, a small, but finite, aver- 



age current I av ~ Ityp/g is obtained [271 ]. Including this 



effect in the framework of this article remains an open 
question. 

Conclusion.- The persistent current has been mostly 
characterized by its first two cumulants. Here, we pro- 
posed a replica trick allowing to calculate at once all 



the cumulants or, equivalcntly, the complete distribution 
function. We mostly applied this trick to the diffusive 
regime, when the statistics is Gaussian and higher or- 
der cumulants are negligible. We believe that the trick 
could be extended to regimes where the Gaussian statis- 
tics breaks down. 

The replica trick introduced in this paper can be ap- 
plied to the evaluation of the probability distribution of 
other thermodynamic quantities. For instance, the non- 
linear sigma-model was used to calculate the mesoscopic 
fluctuations of the supercurrent in metallic Joscphson 
junctions [29[ . We would easily find that the statistics of 
the supercurrent is also Gaussian in the diffusive regime. 
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